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Electron tunneling through a two stage Kondo system constituted by a double quantum-dot 
molecule side coupled to a quantum wire, under the effect of a finite external potential is studied. 
We found that I-V characteristic shows a negative differential conductance region induced by the 
electronic correlation. This phenomenon is a consequence of the properties of the two stage Kondo 
regime under the effect of an external applied potential that takes the system out of equilibrium. 
The problem is solved using the mean- field finite- U slave-boson formalism. 
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I. INTRODUCTION 

Many devices exhibit negative differential conductance 
(NDC) as multiple quantum wells, double barrier, dou- 
ble quantum dots, etoi*2£££&. The NDC has applica- 
tions as amplifiers and oscillators in the microwave, mm- 
wave and Terahertz frequency range. Extensive experi- 
mental and theoretical investigation have been devoted 
to the study of the I-V characteristics and NDC phe- 
nomenon in double quantum dot molecule o 5 ' 6 ' 7 ' 8 '^. More- 
over, there is a wide literature on transport through dou- 
ble quantum dots (DQD) in different geometries, e.g, 
DQD in series and in paralle&i£. Two aspects of elec- 
tronic transport through quantum dots have attracted 
great attention in the last years, the Coulomb block- 
ade effect and the Kondo effect. Recently Kondo effect 
has been studied in side attached 1 ^ and parallel quantum 
dot a 14 i 15 . Recent electron transport experiments showed 
that Kondo and Fano resonances occur simultaneously^. 
Multiple scattering of traveling electronic waves on a lo- 
calized magnetic state are crucial for the formation of 
both resonances. The condition for the Fano resonance 
is the existence of two scattering channels: a discrete 
level and a broad continuum bandit. 

An alternative configuration consists of a double quan- 
tum dot molecule side coupled to a perfect quantum wire 
(QW) 18 . This structure is reminiscent of the T-shaped 
quantum wave guides^. In this case, the QDs act as 
scattering centers in close analogy with the traditional 
Kondo effect. 

Although the electron-electron interaction does play 
an important role in many systems that exhibit nega- 
tive differential conductance, as for instance producing 
bi-stabilities in the current, it is not the driving force of 
the negative conductance itself. In this work instead, we 
study a system in the Kondo regime with a I-V char- 
acteristics that possesses a NDC region that is induced 
by the electronic correlation itself, tuned by the external 
potential and the parameters of the system. 

The system is constituted by a quantum dot molecule 
side coupled to a quantum wire, as shown in Fig. [1] We 
use the finite-C/ slave boson mean-field approach, which 
was initially developed by Kotliar and Ruckenstein^ and 



used later by Bing Dong and X. L. Lei to study the trans- 
port through coupled double quantum dots connected to 
leads^. This approach enforces the correspondence be- 
tween the impurity fermions and the auxiliary bosons to 
a mean-field level to release the U — oo restriction. This 
allows to treat non-perturbatively the dot-lead coupling 
for an arbitrary strength of the Coulomb interaction U 23 . 

In a previous work we study this system in a ther- 
modynamic equilibrium situation. We found that in the 
weak interaction regime, when the direct antiferromag- 
netic interaction between the dots is less than the Kondo 
temperature associated to the internal dot, the transmis- 
sion spectrum shows a structure with two anti-resonances 
localized at the renormalized molecular energies of the 
double quantum do$^. The LDOS of the system shows 
that when the Kondo correlations are dominant the sys- 
tem is in a two stage Kondo regime with two different 
temperatures Tfci,Tfc2, each one associated to a dot. 

In the present paper we study this system under the 
effect of a finite external field, which takes it out of 
the thermodynamic equilibrium, modifying the Kondo 
regime and even destroying it, for enough large fields. 
This process has fundamental consequences on the trans- 
port properties of the system, and in particular, creates a 
remarkable NDC in the I-V characteristics. This NDC 
phenomenon can be understood realizing that the ap- 
plied bias destroys the lower temperature Kondo state of 
the external dot modifying the capability of the other to 
interfere on the current that goes along the system. 



II. MODEL 

Let us consider a quantum dot molecule side coupled to 
a perfect quantum wire (QW) (see Fig. [TJ. We describe 
it by the two-impurity Anderson Hamiltonian. Each dot 
has a single level energy ei (with I = 1,2) and equal intra- 
dot Coulomb repulsion U. The side attached quantum- 
dot molecule is coupled to the QW with coupling to. The 
QW sites have zero local energies and a hopping param- 
eter t. 

The corresponding model Hamiltonian is, 
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FIG. 1: Scheme of double quantum dot (DQD) attached to 
a lead (perfect quantum wire (QW)). The QW is coupled to 
the left (L) and right (R). 
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where c| CT (cj )CT ) is the creation (annihilation) operator of 
an electron with spin a at the i-th site of the quantum 
wire; /, (fl,<r) is the creation (annihilation) operator 
of an electron with spin a in the 1-th QD, ni, a is the 
corresponding number operator, t c is the hopping matrix 
element between the dots and ei correspond to the energy 
of the local states at the dots. 



III. SLAVE BOSON MEAN FIELD THEORY 

To find the solution of this correlated fermions system 
for finite U, we appeal to an analytical approach where, 



generalizing the infinite-t/ slave-boson approximation 21 
the Hilbert space is enlarged at each site, to contain in 
addition to the original fermions a set of four bosons^ 
represented by the creation (annihilation) operators ej 
(ei), p\ a (pi,a), and d\ (di) for the l-th dot. They act as 
projectors onto empty, single occupied (with spin up and 
down) and doubly occupied electron states, respectively. 
Then, each creation (annihilation) operator of an electron 
with spin a in the l-th QD, is substituted by ]\ a Z\ a 
(Zl,afl,a) where: 
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As the problem is solved adopting the U-finite slave bo- 
son mean field approximation (SBMFA), the operator is 
chosen to reproduce the correct C/ — ^ limit in the mean- 
field approximation without changing neither the eigen- 
values nor the eigenvector^ 



The constraint, 
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corporated with Lagrange multipliers A ; and A ; a into 
the Hamiltonian. Also in the mean-field approximation 
all the boson operators are replaced by their expectation 
values pi. ai ei and di which can be chosen, without loss 
of generality, as real numbers. 



The Hamiltonian in this new and enlarged Hilbert 
space, is, H = Hf, + H e , where 
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depends explicitly only upon the boson expectation val- 
ues (e/ = (e;) = ^ e |^ an d equivalently for the others 

operators) and the Lagrange multipliers. The Hamilto- 
nian H e can be written, 
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The tight-binding Hamiltonian depends implicitly on the 
boson expectation values through the parameters: ii_ a — 

(2) ~ ~ ~ ~ 

e i,a + \ t(T i to,a = toZi, a ,tc,a = ic^i,cr^2,ui where Z[^1S 
the value assumed by the operator Zi >a - when the four 
boson operators are substituted by their mean values in 
Equation O 



Pi.cr (ej + di) 



(5) 



The boson operator expectation values and the La- 
grange multipliers are determined by minimizing the en- 
ergy (H) with respect to these quantities. It is obtained 
in this way a set of nonlinear equations for each quantum 
dot, relating the expectation values of the four bosonic 
operators, the three Lagrange multipliers and the elec- 
tronic expectation values, 
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where in the absence of external magnetic field the solu- 
tions are spin independent. 

To obtain the electronic expectation values (•••), the 
Hamiltonian, H e is diagonalized. Their stationary states 
can be written as 
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where a*- and bf are the probability amplitudes to find the 
electron at the site j and at the Z-th QD respectively, with 
energy u = — 2t cos k. As we study the paramagnetic case 
the spin index is neglected. 

The amplitudes a* and bf obey the following linear 
difference equations 
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In order to study the solutions of Eqs. ©, we assume 
that the electrons are described by an incident, a reflected 
and a transmitted plane waves with unitary, r and r am- 
plitudes, respectively. 26 That is, 



e*j + re- ife - J ; (k ■ j < 0) , (9a) 



(fe-j>0). (9b) 



Inserting Eqs. ([9]) into Eqs. ([8]), we get an inhomoge- 
neous system of linear equations for r, r, and 6f, lead- 
ing to the following expression in equilibrium (k = k') 



(w-i-)(u- e+) + i(w - e d2 )T 



(10) 



where the bonding (e_) and antibonding energy (e+) are 
defined by^e± = (e dl + e d2 )/2 ± y/(i dl + e d2 )/2) 2 + l\ 
and r = Trtypo is the renormalized coupling between the 
double quantum-dot and the quantum wire and po is the 
density of states of the leads at the Fermi level. In spite of 
the apparent simplicity of the expression, it is necessary 
to remember that the quantities to and t c implicitly de- 
pend on the expectation values of the boson and fermion 
operators. 

The transmission probability is given by T — |r| 2 , 
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From the amplitudes b\ and b\ we obtain the local 
density of states (LDOS) at the quantum dot I (with 
i = 1,2). 
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In the nonequilibrium case, we suppose a finite source- 
drain biased with a symmetric voltage drop. The incident 
electrons from the left side (L) are in equilibrium with 
the thermodynamical potential pi = P + V/2 and those 
from the right side (R) with pn = p — V/2, 

Once the amplitudes a and bj g are known, the elec- 
tronic expectation values are obtained from, 
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= 2 E ^E/(^-^)^* a h ( 14 ) 

a=L,R k a 

where efc Q = — 2icosfc a . The current given by, 

^^E^^"^) 1 " 1 ^ ^} (15) 

where / (efc a — Ma) it is the Fermi function for incident 
electrons from the a side and the sum on k a is taken up 
to the maximum value cos (— fx a /2t). 
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FIG. 2: LDOS at each dot, left panel QD1, right panel QD2, 
for t c = 0.5V (solid line), V g = — 2F for various values of on 
site energy, U = 6r (solid line), U = 2T (dotted line) U = 
(dashed line) 

The quantities ej, the energies of the local states at the 
dots, are taken to be equal to e\ — e-i = V g , where V g is 
the gate voltage applied to the quantum dots. 

IV. RESULTS 

A. Thermodynamics Equilibrium Case 

In order to obtain a more clear insight, we study first 
the system in thermodynamical equilibrium. The DOS 
at each dot of the quantum-molecule is shown in Fig- 
ure [2 for various values of U, for the case t c = 0.5 T 
and V g = — 2T. As U increases, U > 2T, the system 
passes from the intermediate valence regime into a Kondo 
regime. This process is clearly seen to occur for QD2 
where the resonance shifts and becomes sharper as the 
system enters into the Kondo regime. The same pro- 
cess takes place for the QD1 although the sharpening is 
less accentuated. It is noticeable that as QD2 develops 
a Kondo resonance, a dip appears at the Fermi energy 
in the LDOS of QD1. This is produced because the spin 
of the local electron at QD2 is Kondo correlated with 
the conduction band spins through the mediation of the 
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FIG. 3: LDOS at each dot, left panel QD1, right panel QD2, 
for t c = 0.5F (solid line), U — 4F for various values of the 
gate voltage, V g — — 2r(solid line),F 9 = — T(dotted line) and 
V g = 0(dashed line) 



intermediate dot, that is as well at the Kondo regime. 
This coupling creates a sharp peak at the LDOS of QD2 
with width T^2 and a depletion of this same width at the 
Fermi level of the Kondo peak of width Tui correspond- 
ing to QD1. This result implies the existence of a two 
stage Kondo regime that appears in this system in the 
weak interacting limit, when t 2 c /U < Tki— , being Tki 
and Tfc 2 the Kondo temperatures associated to each dot. 
It is the depletion of the LDOS of QD1 at the Fermi level 
that permits the transmission to be unity. In this case, 
although the side attached QD1, is in the Kondo regime, 
it does not provide an alternative path for the conducting 
electrons. As a consequence, there is no destructive in- 
terference because there is only the direct path available 
for them. 

In Figure [3] we show the LDOS of each QD for different 
values of V g in the two stage Kondo regime. As expected, 
the peaks are pinned at the Fermi energy independent of 
the values of the V g , a clear signal of the Kondo regime. 
Here it is clearly seen at the Fermi energy, the antireso- 
nance at it QD2 and the concomitant resonance at the 
QD1. 



B. The Out of Thermodynamical Equilibrium 
Situation 

The previous discussion was restricted to the thermo- 
dynamical equilibrium situation and was presented with 
the purpose of clarifying the concepts involved. When an 
external potential is applied the scenario changes com- 
pletely. The physics of this new situation can be ex- 
plained analyzing the current and the differential conduc- 
tance dl/dV , two significant and experimentally measur- 
able quantities, as a function of the applied field. 

Figure Q] displays the I-V characteristics (solid line) 
for two values of t c and different values of U, with 
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FIG. 4: 7-V characteristics for t c — 0.5F (solid line) and t c — 
F (dashed line) for various values of U with ei = E2 = —U/2 



V g = —U/2. For [7 = 0, the J-V characteristics shows a 
plateau that is related to the Fano antiresonances in the 
transmission spectrum. When the applied potential is of 
the order of the interdot interaction, the transmission, 
due to the Fano destructive interference, is almost zero 
in the bonding and antibonding regions giving no addi- 
tional contribution to the total current as the applied 
potential is increased. This is the origin of the plateau 
behavior, shown in the Figure EJ when U=0. As U is 
increased, a negative differential conductance appears in 
the I-V characteristics that gets more important as t c is 
augmented. 




FIG. 5: Differential conductance for t c — 0.5F (solid line) and 
t c = r (dashed line), for various values of U with V g = — U/2 

Figure \5\ depicts the differential conductance for the 
same parameters of Figj4] For U = it reflects essen- 
tially the transmission spectrum. As U is increased the 
differential conductance becomes negative in a region of 
the applied potential, reflecting the fact that the NDC 
is a consequence of the Coulomb interaction. It is a small 
effect in the fluctuating valence regime T > U and devel- 
ops completely in the Kondo regime T < U, increasing 
with t r . 



In order to get insight into these results it is convenient 
to write the transmission (Eq.ll) as the superposition of 
a Fano and a Breit-Wigner line shapes, a good approxi- 
mate expression for large values of U— . The results is, 
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where e = ui/T and A — t^/T corresponds to the two 

Kondo temperatures T^i = T and T^2 — A. 

An analytical expression can be obtain for the current 
by integrating over u> the transmission probability given 
in Eq. (H]). 
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We identify each term of the Eq. (fTT)) as follows. The 
first term in the right side is the contribution arising from 
an ideal one dimensional conductor. The second term 
arises from the Kondo-Fano state with temperature Tki 
giving rise to a quasi plateau for the current and almost 
zero differential conductance when \V\ <C T. The third 
term arises from the Kondo state weakly coupled to the 
wire and it is responsible for the rapid increase of the 
current in the region of small applied potentials. 

It is important to emphasize that in this expression the 
quantities A and T are functions of V obtained through 
the self-consistent calculation presented above. The case 
of A = Tk2 is shown in the inset of Fig[51 It is clear 
that increasing V, both the Kondo temperature of the 
external dot and the current reduce to zero simultane- 
ously with the disappearance of the NDC phenomenon. 
As mentioned above, this behavior can be understood 
realizing that the effect of the external quantum dot is 
to reduce the intermediate dot interference effect on the 
current circulating along the leads. This role is exercised 
by the external quantum dot as far as it is at the Kondo 
regime. When the external potential is large enough as 
to disrupt its Kondo ground state (Tk2 < T^)the inter- 
ference is reestablished and the current goes to zero. As 
Tfc2 increases with t c , for larger values of t c this disrupt- 
ing process requires bigger values of the applied poten- 
tial as shown in FigJU When |V| increases still further 
it destroys as well the Kondo regime of the intermediate 
quantum dot and the current rises as depicted in this 
same figure. This seems to be the behavior of a two 
stage Kondo system under the effect of an external po- 
tential. The Kondo regime of the outside quantum dot, 
not directly connected to the continuum, depends upon 
the Kondo effect of the intermediate quantum dot. As 
its Kondo temperature is lower, its Kondo ground state 
is disrupted by lower values of V than the other dot. 
This process manifests in the transport properties by the 
appearance of a NDC region in the current. The FigO 
displays a comparison between the I-V characteristics 
of the numerical calculation and the approximation (Eq. 



(17)) for U = 2T and t c = 0.5T. The approximation over it maintains the form. 

estimates the peak of the current however, qualitatively, Deriving the current in Eq. (|17p we obtain 
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FIG. 6: Comparison of I—V characteristics between the 
numerical calculation (solid line)and the approximation 
(Eq.4)(dashed line) for U = 2r, t c = 0.5T with V g = -U/2. 
T2K = A in the inset 

In the above equation, the last term is responsible for 
the negative differential conductance. This term is neg- 
ative because dA/dV < 0. This expression is able to 
reproduce very well the results of the differential conduc- 
tance shown in FigJS] 

Regarding the observability of the NDC in the side- 
coupled double quantum dot molecule, we consider the 
value of r given by Sato et al^, T = 3meV. For 
tc = 0.5r and U = T the maximum and minimum of the 
current are I m ax ~ 24 nA and I m in ~ 6 nA, respectively, 
giving a peak to valley ratio of 4:1. The lower Kondo 
temperature of the external dot for these same parame- 
ters is of the order of 2 K. These values of the current 



and temperature are well above the experimental limits 
of present day techniques. 



V. SUMMARY 

In summary we have studied the nonequilibrium trans- 
port through a double quantum dot molecule side- 
coupled to a quantum wire using the finite-?/ slave boson 
mean field approach at T — as a function of the pa- 
rameters that define the system. We find that the I-V 
characteristics shows a remarkable NDC, different from 
case reported in the literature, is induced by the elec- 
tronic correlation. This NDC behavior is a consequence 
of the properties of a two stage Kondo system under the 
effect of an external applied potential that takes the sys- 
tem out of equilibrium when the applied potential is large 
enough to destroy the Kondo regime characterized by the 
lower Kondo temperature. 
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